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SUMMARY 


Interfero grams  were  taken  of  the  flow  around  a  sphere  at  Mach  num¬ 
bers  of  1.30  and  1.62.  The  results  of  the  evaluation  of  these  inter - 
ferograms  are  shown  as  plots  of  contours  of  constant  density  ratio  in 
the  flow  field  around  the  sphere  and  as  the  distribution  of  flow  vari¬ 
ables  along  the  axis  of  symmetry  between  the  shock  wave  and  the  sphere. 
In  appendixes  the  theory  of  the  analysis  is  reviewed,  sources  of  error 
are  discussed,  and  the  detailed  procedure  used  in  the  evaluation  of  the 
interferograms  is  described. 


INTRODUCTION 


One  of  the  flow  problems  currently  receiving  attention  is  that  of 
the  flow  behind  detached  shock  waves.  The  present  paper  gives  experi¬ 
mental  data  on  the  flow  behind  the  detached  shock  waves  on  a  sphere. 
These  data  consist  principally  of  density  distributions  around  spheres 
and  were  obtained  for  two  supersonic  Mach  numbers  by  analysis  of  inter¬ 
ferograms.  Previous  investigations  of  axially  symmetric  flow  by  inter¬ 
ferometry  can-- be  found  in  reference  1. 

For  the  convenience  of  those  interested  in  the  application  of 
interferometry  to  the  study  of  axially  symmetric  flow  fields,  the 
symbols  used  are  given  in  appendix  A  and  the  theory  of  the  analysis  is 
reviewed  in  appendix  B.  A  discussion  is  given  in  appendix  C  of  the 
sources  of  error  and  the  accuracy  of  the  results  obtainable  from  inter¬ 
ferograms  of  axially  symmetric  flow.  (This  discussion  is  somewhat 
similar  to  those  given  in  references  1  and  2.)  Because  reasonable 
accuracy  in  the  results  requires  that  much  care  be  taken  in  measuring 
and  in  plotting  fringe  shifts,  a  detailed  description  of  the  computa¬ 
tional  procedure  that  was  used  to  obtain  the  final  results  is  given  in 
appendix  D.  A  table  of  the  coefficients  for  5 0  zones  for  use  in  the 
Weyl  method  of  evaluation  has  been  calculated  and  is  also  given. 
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APPARATUS 


The  test  setup  consisted  principally  of  a  supersonic- jet  apparatus 
that  gave  a  3-inch-square  free  jet  and  an  interferometer  with  4-inch- 
square  plates.  This  apparatus  has  been  described  in  reference  3.  The 
only  change  in  the  apparatus  was  a  change  in  the  film  magazine  to  permit 
the  use  of  70-millimeter  film  rather  than  the  35-millimeter  size  previ¬ 
ously  used.  Two  different  Mach  numbers  were  obtained  by  using  two 
supersonic  nozzles  similar  to  those  described  in  reference  3»  The  test 
model  was  a  steel  sphere  l/2  inch  in  diameter  that  was  supported  from 
behind  by  a  sting  l/8  inch  in  diameter. 


RESULTS  AND  DISCUSSION 
Interferograms 


An  interferogram  that  showed  the  undisturbed  fringes  in  the  test 
section  was  taken  with  no  flow  (fig.  l)  and  then  an  interferogram  of 
the  flow  around  the  sphere  at  Mach  number  1.30  was  taken  (fig.  2).  The 
same  procedure  was  followed  at  a  flow  Mach  number  of  1.62}  the  no-flow 
interferogram  is  shown  in  figure  3  and  the  flow  interferogram,  in  fig¬ 
ure  4.  Cross  sections  were  drawn  on  the  flow  interferograms  perpendicu¬ 
lar  to  the  axis  of  symmetry.  The  variation  of  density  along  each  cross 
section  was  then  obtained  by  the  method  described  in  appendixes  B  and  D. 
A  plot  of  the  ratio  of  density  along  the  cross  section  to  free-stream 
density  was  made  for  each  cross  section. 


Density  Distributions 

The  plots  of  density  ratio  along  the  cross  sections  were  used  for 
plotting  contours  of  constant  density  ratio  in  the  flow  field  around 
the  sphere.  Experimental  points  taken  from  the  unfaired  plots  of  den¬ 
sity  ratio  are  shown  in  figure  5*  The  curves  that  are  shown  are  faired 
through  these  points. 

The  method  of  analysis  of  the  interferograms  that  was  used  has  the 
disadvantage  that  the  results  obtained  are  inaccurate  in  a  small  region 
immediately  behind  a  discontinuity  in  the  density.  (The  reason  for  the 
inaccuracy  is  explained  in  appendix  C.)  Because  of  this  inaccuracy  in 
the  results  immediately  behind  the  shock  wave,  points  on  the  plots  of 
density  ratio  that  were  near  the  shock  wave  and  that  were  obviously  in 
error  were  not  included  in  figure  5*  For  each  contour  in  figure  5, 
therefore,  a  small  region  ^adjacent  to  the  shock  wave  exists  in  which 
experimental  points  have  not  been  obtained.  The  density  contours  have 
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not  been  extended  through  this  region  to  the  shock  wave .  On  the  shock 
wave,  however,  the  correct  location  of  the  intersection  of  each  contour 
with  the  shock  wave  has  been  indicated.  The  locations  of  the  inter¬ 
sections  were  obtained,  from  the  equation  for  the  density  ratio  across  a 
two-dimensional  oblique  shock  wave  with  use  of  measured  slopes  of  the 
shock  wave. 

The  magnitude  of  any  systematic  errors  in  the  contours  shown  in 
figure  5  is  not  known.  The  uncertainty  in  the  location  of  the  contours, 
aside  from  any  systematic  errors,  is  believed  to  be  not  larger  than  plus 
or  minus  approximately  one-half  the  distance  between  adjacent  contours. 


Flow  Variables  along  Axis 

The  distribution  of  the  flow  variables  along  the  axis  of  symmetry 
is  of  interest.  The  region  of  the  axis  is,  however,  a  region  of  rela¬ 
tively  large  error.  In  order  to  obtain  what  are  thought  to  be  the  best 
values  of  density  ratio  along  the  axis,  the  values  at  the  axis,  obtained 
by. fairing  the  plots  of  density  ratio  for  each  cross  section,  were  used 
to  plot  the  variation  of  density  ratio  along  the  axis  between  the  shock 
wave  and  the  sphere  (fig.  6).  Immediately  behind  the  normal  part  of  the 
shock  wave,  the  density  ratio  is  well-known  from  the  Rankine -Hugoniot 
relations.  At  the  nose  of  the  sphere,  the  correct  value  of  the  stagna¬ 
tion  density  is  also  known.  These  theoretical  values  were  used  in 
fairing  the  curves  shown  in  figure  6  through  the  experimental  points. 

The  density  ratios  along  the  axis  given  by  the  faired  curves  have  also 
been  converted  to  ratio  of  static  pressure  to  stagnation  pressure 
(fig.  7).  The  distribution  of  density  ratio  along  the  axis  has  also 
been  converted  into  the  distribution  of  Mach  number  shown  in  figure  8. 
(The  stagnation  temperature  for  the  free-stream  Mach  number  of  1.30 
was  524. V3  F  abs.  and"  for  Mach  number  1.62  was  533.1°  F  abs.  The  free- 
stream  density  for  Mach  number  1.30  was  0.003215  slug/cu  ft  and  for 
Mach  number  1,62  was  0.003501  slug/cu  ft.) 

The  variation  of  density  ratio  along  the  surface  of  the  sphere  for 
Mach  number  1.62  is  shown  in  figure  9*  The  experimental  points  were 
taken  from  the  plots  of  density  ratio  along  each  cross  section. 


CONCLUDING  REMARKS 


Experimental  data  on  the  flow  around  spheres  at  supersonic  speeds 
have  been  presented.  Although  greater  accuracy  in  the  data  would  be 
desirable,  the  data  given  herein  are  accurate  enough  that  they  should 
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be  of  use  in  checking  and  evaluating  theoretical  treatments  of  the 
problem  of  predicting  the  flow  behind  detached  shock  waves. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base,  Va. ,  June  12,  1950 
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APPENDIX  A 

'SYMBOLS 

* 

radius  of  first  zone 
coefficient 

axial  distance  from  center  of  curvature  of  shock  wave 
number  of  zone  at  which  density  is  required 
zone  number 

distance  through  test  section  or  through  disturbance 
Mach  number 

number  of  outermost  zone 
index  of  refraction 
pressure 

radial  distance  from  axis  of  symmetry 
radius 

fringe  shift  per  unit  fringe  width 
temperature 

value  of  v  for  which  density  is  required 

Cartesian  coordinate  in  direction  of  undeviated  light  ray 

Cartesian  coordinate  perpendicular  to  light  ray  and  to  free- 
stream  flow  direction 

zone  width 
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ratio  of  specific  heats 
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k  Gladstone -Dale  constant 

X0  wavelength  of  light  in  vacuum 

V  difference  between  index  of  refraction  at  point  r  and 

reference  index  of  refraction 

p  density 

p  reference  density 

ref 

0  angle  between  ray  and  normal  to  surface 

Subscripts: 
s  stagnation 

0  free  stream 

1  unrefracted  ray;  also,  correct  value 

2  refracted  ray;  also,  approximation  to  correct  value 
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APPENDIX  B 

ANALYSIS 


The  problem  is  to  determine  the  density  from  the  known  fringe 
shift.  For  the  two-dimensional  case,  where  the  density  is  not  a  func¬ 
tion  of  x,  the  relation  between  fringe  shift  and  density  has  been 
shown  in  reference  3;  for  example,  to  be 


For  the  three-dimensional  case,  where  the  density  is  a  function  also 
of  x,  the  fringe  shift  at  any  point  on  the  interferogram  is  the  sum¬ 
mation  of  the  shifts  produced  in  each  element  of  path  of  the  ray  of 
light.  If  the  assumption  is  made  that  the  density  distribution  is 
axially  symmetric,  then  p  is  a  function  of  r  alone.  If  the  addi¬ 
tional  assumption  is  made  that,  within  the  disturbance,  the  light  ray 
does  not  deviate  from  a  straight  line  that  is  parallel  to  the  x-axis, 
then  the  relation  between  fringe  shift  and  density  is 


Three  methods  of  handling  equation  (l)  have  been  proposed.  The 
earliest  is  that  which  is  discussed  by  Schardin  (reference  4).  In  this 
method,  the  disturbance  is  divided  into  zones,  as  shown  in  figure  10, 
where  a  zone  is  the  region  between  two  adjacent  concentric  circles. 

(The  zone  between  circles  i  and  i  +  1,  for  example,  is  referred  to 
in  the  present  analysis  as  zone  i.)  The  assumption  is  then  made  that 
the  density  is  constant  in  each  zone.  The  resulting  equation,  which  is 
thus  an  approximation  to  equation  (l),  can  be  solved  for  p(r). 

In  the  second  method,  which  was  introduced  by  Van  Voorhis  (refer¬ 
ence  1),  the  value  of  the  quantity  p(r)  -  Pref  in  each  zone  is  assumed 

to  be  a  linear  function  of  the  radius  r  of  the  zone.  The  resulting 
equation,  which  is  thus  also  an  approximation  to  equation  (l),  is 
solvable . 

In  the  third  method,  which  was  introduced  by  Weyl  (reference  2), 
equation  (l)  is  solved  exactly  for  p(r)  -  Pref  and  'the  resulting 

definite  integral  is  approximated  by  assuming  that  in  each  zone  the 
fringe  shift  is  a  linear  function  of  y2.  First  the  variables 
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and 


are  substituted  into  equation  (l)  to  give  the  equation 


x  1  /'"v  (v)  dv 

u>  •  ~ol  WT 

This  integral  equation  is  the  same  as  Volterra's  integral  equation  of 
the  first  kind  except  that  the  limits  of  integration  differ  from  the 
Volterra  equation  limits  of  0  and  u.  This  equation  can  be  solved 
in  a  manner  similar  to  that  used  for  solving  Volterra1 s  equation  (refer¬ 
ences  2  and  5)*  The  result  is 


•v(w) 


I"'00  S'(u)  du 
71  J  V  \ju  -  w 


where  the  prime  indicates  differentiation  with  respect  to  u.  This 
resulting  integral  can  be  computed  by  breaking  the  integral  into  steps, 
one  step  for  each  zone,  and  by  assuming  that  in  each  zone  the  fringe 

shift  S(u)  is  a  linear  function  of  u,  where  ri2  <u<  ri+]_2*  Then 
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For  zones  of  equal  width  z. 


and 


ri  =  ^-z 


2X0Jtl 


Jtz 


H(« 


k=i 


o  \  V(k  +  l)2  -  i2  -\/k2'-  i2 

Sk+l)  2k  +  1 - 


(2) 


Equation  (2)  is  an  approximation  to  the  solution  of  equation  (l)  and 
is  the  relation  that  was  used  for  computing  the  results  given  in  the 
present  paper. 

The  coefficients 

A  _  \/(k  +  1  )2  ZZ  -  V^2  -  i2 

k,i  2k  +  1 

for  50  zones  were  calculated  by  the  Bell  Telephone  Laboratories 
X-66jkk  relay  computer  at  the  Langley  Laboratory  and  are  given  in 
table  I. 

The  situation  may  arise  that  would  make  desirable  the  use  of  more 
than  50  zones.  Because  the 'number  of  coefficients  is  approximately 
equal  to  half  the  square  of  the  number  of  zones,  a  large  number  of  coeiV 
ficients  would  have  to  be  calculated.  A  method  has  been  devised  in 
which  the  number  of  coefficients  is  equal  only  to  the  number  of  zones 
and  each  coefficient  is  also  easier  to  calculate.  The  idea  is  to  use 
zones  not  of  equal  width  but  of  radius  proportional  to  the  square  root 
of  i, 

r  =  a\fT 

z  =  ri+l  "  ri  * 

-  a(\fi  +  1  -  \|i) 

where  a  is  the  radius  of  the  first  zone.  Then 


v  -?  = 


Tta 


y'1  (sk  -  sk+l)  ( '\fk~-  i  +  1  -  \/k  -  i) 


k=i 
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The  coefficient  is  a  function  only  of  k  -  i,  which  ranges  in  value 
from  0  to  N  -  1.  Therefore,  only  N  coefficients  need  be  calculated, 
rather  than  N2/2.  This  new  method  was  not  used  in  the  calculations  for 
the  present  paper  but  is  mentioned  because  of  its  possible  value  in 
future  calculations. 
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APPENDIX  C 

SOURCES  OF  ERROR 


The  purpose  of  the  present  section  is  to  consider  some  of  the 
sources  of  error  inherent  in  the  application  of  equation  (2).'  Equa-  ■ 
tion  (2)  can  be  used  for  determining  density  distribution  with  suffi¬ 
cient  accuracy  only  if  the  following  conditions  are  sufficiently  well 
satisfied: 

.  (a)  The  index  of  refraction  of  the  gas  varies  linearly  with  density 

(b)  A  reference  density  at  some  point  in  the  field  of  view  of  the 
interferometer  can  be  determined 

(c)  The  disturbance  is  rotationally  symmetrical  about  an  axis  that 
is  known 

(d)  The  light  beam  is  not  refracted 

(e)  The  variation  with  y  of  fringe  shift  S(y)  for  a  cross 
section  of  the  disturbance  perpendicular  to  the  axis  of  symmetry  is 
known 

(f)  The  variation  with  u  of  fringe  shift  S(u)  can  be  replaced 
by  a  series  of  straight  lines 

As  far  as  is  known,  condition  (a)  is  always  sufficiently  well 
satisfied  by  air.  Whether  condition  (b)  is  satisfied  at  all  is  gener¬ 
ally  determined  by  the  experimental  setup.  For  the  present  investiga¬ 
tion  the  reference  density  was  the  free-stream  density.  For  determining 
the  free-stream  density,  the  value  of  the  stagnation  temperature  was 
measured  at  the  instant  the  interferogram  was  taken.  The  Mach  number 
given  by  the  nozzle  was  known  from  a  previous  calibration  described  in 
reference  6.  The  free-stream  temperature  was  then  calculated  from  the 
relation 


Then  the  free-stream  density  was  computed  by  the  general  gas  law  under 
the  assumption  that  the  free-stream  static  pressure  was  equal  to  the 
ambient  atmospheric  pressure,  which  was  measured.  The  determination  of 
free-stream  density  is,  of  course,  subject  to  some  degree  of  uncertainty. 
The  error  in  free-stream  density  is  estimated  to  be  not  greater  than 
3  percent,  and  the  probable  error  is  estimated  to  be  about  1  percent. 
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Condition  (c)  is  assumed  to  t>e  satisfied  for  the  part  of  the  flow 
field  that  is  behind  the  shock  wave  and  ahead  of  the  most  upstream  part 
of  the  intersection  of  the  shock  wave  with  the  mixing  region  at  the 
edges  of  the  jet.  Because  the  shock  wave  is  axially  symmetric  and  the 
jet  is  square,  as  soon  as  the  shock  wave  reaches  the  edge  of  the  jet, 
condition  (c)  is  no  longer  satisfied. 

Condition  (d),  that  the  light  beam  is  not  refracted,  is  usually 
sufficiently  well  satisfied  in  rotationally  symmetric  flow.  (The  only 
place  the  effects  of  refraction  are  likely  to  be  large  is  in  boundary 
layers  in  two-dimensional  flow.)  That  the  effect  of  refraction  on 
fringe  shift  and  on  calculated  density  was  fairly  small  in  the  present 
case  was  determined  by  calculating  the  effect.  The  purpose  of  the 
calculation  was  not  to  make  a  precise  analysis  of  the  effects  of  refrac¬ 
tion  but  to  determine  approximately  the  magnitude  of  the  effect  on  cal¬ 
culated  density  in  the  region  where  the  effect  was  largest.  That  region 
is  at  the  nose  of  the  shock  wave.  Certain  assumptions  that  greatly 
simplified  the  analysis  were  therefore  justified.  The  assumption  was 
made  that  the  region  behind  the  nose  of  the  shock  wave  was  a  spherical 
region  of  constant  density.  When  a  ray  of  light  impinges  obliquely  on 
a  region  of  greater  density,  the  ray  is  refracted  and  emerges  at  a 
place  different  from  the  place  it  would  emerge  if  not  refracted.  The 
optical  path  length  is,  therefore,  changed  by  the  refraction.  Because 
of  the  change  in  optical  path,  the  fringe  shift  produced  by  the  ray  is 
also  changed.  Figure  11  shows  the  ray  of  light  entering  the  sphere  . 
at  Pi  and  leaving  at  P2.  Only  a  meridian  plane,  defined  by  the  ray 
and  the  center  of  the  sphere,  need  be  considered.  The  geometrical  path 
length  without  refraction  is  Lp  and  with  refraction  is  Lg 

Lp  =  2R  cos  0p 

=  2R\]l  -  sin^0p 

Lg  =  2R  cos  02 
=  2R\jl  -  sin^02 

sin  01  =  | 

and,  inasmuch  as 


n2  sin  02  =  ^  sin  0^ 
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then 

nld 

2  ”  r^R 

Therefore, 


The  value  of  np  outside  the  sphere  (in  the  free  stream)  is  1.0004  at 
Mach  number  1.62,  The  density  within  the  sphere  is  assumed  to  be  twice 
that  outside  the  sphere.  Then,  because  n  -  1  is  proportional  to  the 
density,  the  value  of  n2  inside  the  sphere  is  1.0008.  The  ratio  of  n2 
to  np  is  1.0004.  The  values  of  L2/Lp  for  various  values  of.  d/R 
are  shown  in  the  following  table,  together  with  the  percentage  error  in 
the  fringe  shift.  The  corresponding  values  of  distance  behind  the  nose 
of  the  assumed  sphere  of  constant  density,  expressed  in  terms  of  the 
distance  between  the  nose  of  the  shock  wave  and  the  nose  of  the  sphere, 
are  shown  in  the  last  column  of  the  table. 


sin  0 


d/R 

— 

1 — 1 

cvT 

Error  in 
fringe  shift, 
percent 

Percent  of 
distance  from  shock 
wave  to  sphere 

0.995 

1.039 

3.9 

3 

.990 

1.019 

'1.9 

5 

.980 

1.010 

1.0 

11 

.960 

1.005 

•  5 

22 

.9^0 

1.003 

.3 

32  • 

.920 

1.002 

.2 

^3 

The  values  of  the  radius  of  curvature  of  the  shock  wave  at  its  nose, 

2.9  times  the  radius  of  the  model,  and  of  the  axial  distance  between 
the  shock  wave  and  the  model,  0,56  times  the  radius  of  the  model,  were 
taken  from  reference  6. 

The  effect  of  refraction  is  seen  to  be  small  except  immediately 
behind  the  shock  wave.  There,  however,  the  error  due  to  refraction  may 
still  be  relatively  unimportant  because  at  the  same  place  there  may  be 
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a  larger  error  due  to  the  "over shooting”  in  the  computed  density  that 
occurs  immediately  behind  a  discontinuity  in  the  density,  such  as  exists 
at  a  shock  wave.  (The  overshooting  is  explained  subsequently  in  the 
discussion  of  condition  (f).) 

Whether  condition  (e)  is  sufficiently  well  satisfied  for  the  com¬ 
puted  value  of  density  to  lie  within  given  limits  of  uncertainty  depends 
on  a  number  of  factors.  A  relation  that  gives  the  effect  of  some  of 
these  factors  is  derived.  Figure  10  shows  a  cross  section  of  the  flow 
field  divided  into  concentric  zones.  A  typical  ray  of  light  is  shown 
traversing  the  field.  The  innermost  zone  through  which  this  ray  travels 
i$  designated  zone  i.  The  assumption  is  then  made  that,  for  this  ray, 
the  fringe  shift  due  to  the  density  outside  of  zone  i  is  accurately 
known.  The  problem  is  to  determine  how  much  error  in  the  fringe  shift 
at  y  =  i  can  be  tolerated  for  a  given  uncertainty  in  the  determination 
of  the  density  at  the  same  place.  Let  S be  the  correct  value  of 

fringe  shift  at  y  =  i  and  the  approximation  to  •  The  error 

in  fringe  shift  is 


AS  =  Sii  -  Si2 

From  the  definition  of  v  and  equation  (2), 

vix  =  K(Pl  -  pref) 


Then,  because  all  the  fringe  shifts  from  k  =  i  +  1  to  k  =  N  -  1  are 
assumed  to  be  accurately  known  and  are  the  same  for  both  case  1  and 
case  2, 


Ap  =  px  -  p2 


JCZk 


AS  A, 


1*1 


(3) 
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Equation  (3)  can  tie  used  to  calculate  the  accuracy  with  which  fringe 
shift  must  be  known  for  determining  density  within  specified  limits. 

The  equation  is  used  herein  to  investigate  the  effect  of  error  in 
the  value  of  the  fringe  shift  on  the  value  of  the  density  for  a  sample 
case.  A  square  density  distribution  as  shown  in  figure  12  is  assumed. 
The  corresponding  fringe-shift  curve  is  also  shown.  The  assumption  is 
made  that  this  fringe-shift  curve  is  known  to  ±0.001.  From  this  fringe- 
shift  curve  the  density  distribution  is  in  turn  recalculated.  The 
result  is  the  solid  curve  that  peaks  at  zone  24.  The  overshooting  that 
occurs  at  the  discontinuity  in  density  is  not  the  result  of  the  uncer¬ 
tainty  of  0.001  in  the  fringe  shift  at  each  zone  but  is  caused  by  fail¬ 
ure  of  condition  (f)  to  be  satisfied  and  is  discussed  subsequently. 

For  comparison,  the  assumptions  are  also  made  that  the  fringe  shift 
is  0.02  too  high  at  every  zone  and  also  that  the  fringe  shift  is  0.02 
too  low  at  every  zone.  The  density  distributions  calculated  under  these 
assumptions  are  not  shown  in  figure  12  since  they  are  almost  indistin¬ 
guishable  from  each  other  and  also  from  the  density  distribution 
obtained  under  the  assumption  of  an  uncertainty  of  0.001. 

Equation  (3)  is  now  applied  and  is  found  to  give,  at  zone  15,  for 
AS  =  0.02,  a  value  of  Ap  of  0.00014  slug  per  cubic  foot.  The  actual 
value  of  Ap  is  0.00001  slug  per  cubic  foot.  For  this  case,  therefore, 
equation  (3)  gives  a  very  conservative  estimate  of  the  error  in  density. 

Applying  equation  (3)  to  a  case  where  it  will  give  the  least  con¬ 
servative  estimate  of  error  is  desirable.  The  largest  error  in  density 
occurs  when  the  error  in  the  fringe-shift  curve  is  an  oscillating  one 
that  changes  sign  at  each  zone.  The  density  distribution  obtained 
under  this  assumption  of  fringe-shift  error  is  shown  by  the  oscillating 
curve  in  figure  12  for  |AS|  =  0.02.  At  zone  15  the  error  in  density 
is  0.0001  slug  per  cuMc  foot.  This  value  is  to  he  compared  with  the 
value  0.00014  slug  per  cubic  foot  calculated  by  equation  (3)*  Equa¬ 
tion  (3),  therefore,  appears  to  be  a  good  approximation  in  this  case, 
which  is  one  that  simulates  those  that  occur  in  actual  practice. 

Satisfying  condition  (e),  that  is,  obtaining  fringe  shift  with 
sufficient  accuracy,  is  perhaps  the  most  difficult  requirement  to  meet 
in  obtaining  density  distributions  in  three-dimensional  flows.  Because 
equation  (3)  so  clearly  illustrates  the  difficulty,  the  equation  is 
used  in  another  example  to  show  how  accurately  the  fringe  shifts  must  - 
be  known  for  the  density  to  be  determined  accurately.  A  typical  posi¬ 
tion  in  figure  5(b)  is  chosen  and  the  requirement  is  made  that  the 
density  be  obtained  within  1  percent.  A  position  is  chosen  about  midway 


1 6 


MCA  TN  2173 


along  the  contour  on  which  the  ratio  of  density  to  free-stream  density 
is  2.0.  For  this  position, 

_5 

XQ  =  5.17  X  10  centimeters 
z  =  2.3  X  10"^  centimeters 
k  =  O.II67  cubic  foot  per  slug 
i  =  2b 

Ai,i  =  0-14 

p  =  0.0070  slug  per  cubic  foot 
A p  =  0.00007  slug  per  cubic  foot 
Then,  by  equation  (3); 

AS  =  0.04 


In  the  present  investigation,  fringe  shifts  could  generally  be  deter¬ 
mined  with  approximately  twice  this  amount  of  uncertainty,  or  AS 
approximately  equal  to  0.1,  for  the  interferogram  of  the  flow  at  Mach 
number  1.62.  For  the  interferogram  taken  at  Mach  number  1.30,  howevei 
the  fr inge  shifts  were  somewhat  more  uncertain  than  for  the  other  inte: 
ferogram.  The  model  is  farther  from  the  end  of  the  nozzle,  and  turbu¬ 
lence  in  the  mixing  regions  on  the  two  sides  of  the  jet  through  which 
the  light  passes  causes  appreciable  irregularities  in  the  fringe  patter 

That  the  effect  of  shrinkage  of  the  film  during  development  was 
negligibly  small  was  found  by  photographing  a  transparent  scale. 
Furthermore,  by  enlarging  photographs  of  a  scale,  it  was  found  that  dis¬ 
tortions  introduced  by  the  lens  of  the  enlarger  were  also  negligibly 
small . 

Condition  (f)  is  that  S(u)  can  be  replaced  by  a  straight  line  in 
each  zone.  Weyl  has  shown  (reference  2)  that  this  condition  is  suffi¬ 
ciently  well  satisfied  everywhere  except  in  zones  in  which  the  rate  of 
change  of  Sf(u)  is  large,  that  is,  where  the  second  derivative 
of  S(u)  has  large  absolute  values.  At  shock  waves  the  second  deriva¬ 
tive  of  S(u)  is  very  large.  The  result  is  that  condition  (f)  is 
grossly  violated  and  the  calculated  density  overshoots  at  shock  waves. 

At  other  places  where  the  second  derivative  is  large,  but  not  so  large 
as  at  shock  waves,  condition  (f)  can,  by  the  use  of  smaller  zone  widths, 
generally  be  sufficiently  well  satisfied  in  that  the  error  due  to 
violation  of  condition  (f)  is  acceptably  small. 
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APPENDIX  D 


COMPUTATIONAL  PROCEDURE 

The  negatives  of  the  interfero grains  were  about  inches  in  diame¬ 
ter.  Prints  were  made  in  which  the  negatives  were  enlarged  about 
16  diameters  and  these  prints  were  used  for  all  measurements. 

The  computational  procedure  used  is  as  follows: 

The  axis  of  symmetry  was  drawn  on  each  interferogram.  On  each  no¬ 
flow  interferogram  the  fringes  were  numbered  arbitrarily.  A  cross  sec¬ 
tion  was  drawn  perpendicular  to  the  axis  of  symmetry  and  the  position 
of  each  fringe  was  measured  along  the  cross  section.  Then  fringe  number 
was  plotted  as  a  function  of  position.  A  large  scale  that  will  permit 
reading  the  curve  with  the  required  accuracy  had  to  be  used. 

On  the  flow  interferograms,  lines  perpendicular  to  the  axis  of  sym¬ 
metry  were  drawn.  These  lines  were  spaced  about  l/lO  inch  apart  and 
denoted  the  cross  sections  along  each  of  which  the  fringe  shift  was  to 
be  measured.  Again  the  fringes  were  arbitrarily  numbered  and  the  posi¬ 
tions  of  the  fringes  were  measured  along  each  cross  section.  The  posi¬ 
tions  of  both  the  black  and  the  white  fringes  were  measured,  because  on 
part  of  the  interferogram  the  fringes  run  almost  parallel  to  the  cross 
sections  and  the  number  of  points  is,  therefore,  less  than  is  desirable. 

The  next  step  was*  to  determine  the  fringe  shift  at  each  measured 
point  on  each  cross  section.  The  position  of  each  fringe  is  known.  The 
position  of  the  disturbed  fringe  was  used  in  the  previously  mentioned 
plot  of  fringe  number  as  a  function  of  fringe  position  to  determine  the 
number  of  the  undisturbed  fringe  at  that  position.  Then  the  number  of 
the  undisturbed  fringe  was  subtracted  from  the  number  of  the  disturbed 
fringe  to  give  the  fringe  shift. 

Because  the  fringes  had  been  arbitrarily  numbered,  all  fringe 
shifts  were  adjusted  by  the  amount  necessary  to  make  the  undisturbed 
fringe  at  the  same  location  in  both  flow  and  no-flow  interfero grams 
have  the  same  number.  This  adjustment  is  tantamount  to  numbering  the 
fringes  on  either  the  flow  or  the  no-flow  interferogram  not  arbitrarily 
but  in  such  a  way  that  the  same  number  is  given  to  the  same  position  in 
the  undisturbed  portion  of  both  interf erograms . 

Because  the  free-stream  density  was  to  be  used  as  the  reference 
density  in  computing  density  distribution  around  the  sphere  and  because 
the  fringe  shift  must  be  assumed  to  be' zero  for  the  reference  density. 
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the  fringe  shift  in  the  free  stream  had  to  be  determined  and  subtracted 
from  all  other  fringe  shifts.  The  fringe  shift  in  the  free  stream  was 
determined  from  the  free-stream  density  which  was  calculated  as  described 
in  appendix  C. 

A  cross  section  in  the  free  stream  was  chosen  as  giving  the  posi¬ 
tions  of  the  reference  fringes.  The  fringes  in  the  free  stream  were 
at  a  slight  angle  to  the  axis  and,  therefore,  to  the  cross  sections 
which  were  drawn  perpendicular  to  the  axis.  This  tilt  of  the  free- 
stream  fringes  has  an  effect  on  the  fringe  shifts  measured  along  each 
cross  section  in  that  it  causes  them  to  be  different  by  an  amount  that 
is  proportional  to  the  tangent  of  the  angle  of  tilt  and  to  the  distance 
of  the  cross  section  from  the  reference  cross  section.  Proper  adjust¬ 
ment  of  the  measured  fringe  shifts  along  each  cross  section  was  made  to 
eliminate  the  effect  of  the  tilt. 

The  next  step  was  to  plot  fringe  shift  as  a  function  of  position. 

Here  again,  care  must  be  taken  to  use  a  sufficiently  large  scale.  The 
plot  was  then  divided  into  zones,  with  the  zeroeth  zone  at  the  axis  of 
symmetry.  Zone  width  is  generally  determined  by  a  number  of  factors. 
Equation  (3)  shows  that  the  greater  the  zone  width  the  less  the  error 
in  density  due  to  violation  of  condition  (e).  On  the  other  hand,  the 
greater  the  zone  width,  the  smaller  the  number  of  zones,  the  less  well 
condition  (f )  is  satisfied,  and  the  smaller  the  number  of  experimental 
points  that  will  be  obtained.  A  compromise  in  zone  width  is,  therefore, 
necessary.  For  the  present  investigation,  the  cross  sections  were 
separated  into  groups  of  nine  each.  The  same  zone  width  was  used  for 
all  cross  sections  in  a  group  and  was  such  that  the  longest  cross  sec¬ 
tion  in  each  group  was  divided  into  50  zones. 

The  next  step  was  to  read  the  fringe  shift  at  each  zone  and  then 
to  determine,  by  subtracting,  the  change  in  fringe  shift  between  adja¬ 
cent  zones.  This  last  quantity  is  the  quantity  S^.  -  S^+q  that  is 

used  in  equation  (2)  for  calculating  v.  The  computations  were  made  by 
the  Bell  computer.  From  the  values  of  v  the  values  of  the  ratio  of 
the  density  to  the  free-stream  density  were  computed.  These  values  were 
used  to  obtain  a  plot  of  the  density  ratio  against  zone  number  for  each 
cross  section.  A  typical  plot  is  shown  in  figure  13,  which  is  for  a 
cross  section  that  is  about  three-fourths  of  the  way  between  the  nose 
of  the  shock  wave  and  the  nose  of  the  sphere. 

The  plot  for  each  cross  section  has  two  regions  where  the  inac¬ 
curacy  or  uncertainty  is  greater  than  for  the  rest  of  the  curve.  One 
of  these  regions  is  at  and  near  the  shock  wave.  There  the  inaccuracy 
is  caused  by  failure  to  satisfy  condition  (f),  and  the  result  is  over¬ 
shooting  or  undershooting  at  the  shock  wave.  (in  fig.  13  the  shock  wave 
is  at  zone  43  and  the  correct  density  ratio  behind  the  shock  wave  is  1.77*) 
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The  density  ratio  across  a  shock  wave  is  known,  however,  from  well- 
established  theory.  The  correct  values  of  the  density  ratio  behind  the 
shock  wave,  therefore,  can  be  shown  on  the  plots  of  density  ratio.  The 
other  region  of  relatively  large  inaccuracy  is  at  and  near  the  axis  of 
symmetry  (zones  0  to  5  in  fig.  13).  That  the  region  of  the  axis  is  a 
region  of  relatively  large  inaccuracy  can  be  seen  from  equation  (3), 
which  shows  that  the  uncertainty  in  density  is  proportional  to  the  coef¬ 
ficient  A-j^.  The  coefficient  A^  ^  can  be  seen  from  table  I  to  be 

larger  for  zone  0,  which  is  at  the  axis  of  symmetry,  than  for  any  other 
zone.  Unlike  the  other  region  of  relatively  large  inaccuracy,  however, 
for  which  the  density  ratio  immediately  behind  the  shock  wave  is  known, 
the  density-ratio  distribution  along  the  axis  of  symmetry  is  not  known 
from  theory. 
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TABLE  I.-  COEFFICIENTS 


k,  i 


v/(k  +  i)g  -  ±g  -  i/g~T 
2k  +  1 


0  1.00000 
1  0.33333  0.57735 


2 

.20000 

.21928 

0.44721 

3 

.14286 

.14922 

.175^3 

0.37796 

4 

.11111 

.n4oo 

.12427 

.15047 

0.33333 

5 

.090909 

.092464 

.097662 

.10874 

.13383 

0.30151 

6 

.076923 

.077856 

.080873 

•.O868OO 

.097879 

.12172 

0.27735 

7 

.066667 

.067270 

.069184 

.072776 

.078909 

.089735 

.11240 

0.25820 

8 

.058824 

.059236 

.060529 

.062887 

.066709 

.072842 

.083335 

.10493 

0.24254 

9 

.052632 

.052926 

.053842 

.055480 

.058047 

.061944 

.067989 

.078135 

.098784  0.22942 

10 

.047619 

.047837 

.048509 

.049696 

.051514 

.054176 

.058074 

.063993 

.073802 

.093603  0.21822 

n 

.043478 

.043644 

.044153 

.045041 

.046381 

.048294 

.050989 

,05*i848 

. 060628 

,070117 

.089160 

12 

,o4oooo 

.040129 

.040523 

. 041206 

.042224 

.043652 

.045610 

.048306 

.052107 

• 0577^3 

.066935 

13 

.037037 

.037139 

.037451 

.037988 

.038781 

.039878 

.041354 

.043330 

.046007 

.049740 

.055235 

14 

.034^3 

.034565 

.034B16 

.035247 

.035877 

.036739 

.037883 

.039384 

.041360 

.044007 

.047668 

15 

.032258 

.032325 

.032530 

.032880 

.033390 

.034082 

.034989 

.036161 

.037672 

.039637 

.042247 

1 6 

.030303 

.030359 

.030528 

.030817 

.031236 

.031800 

.032532 

.033468 

.034654 

.036165 

.038113 

17 

.028571 

.028618 

.028760 

.029001 

.029349 

0 029815 

.030417 

.03U77 

.032129 

.033321 

.034826 

18 

.027027 

.027067 

.027186 

.027390 

.027682 

.028072 

.028573 

.029200 

.029978 

.030939 

.032131 

19 

.025641 

.025675 

.025777 

.025950 

.026198 

0 026529 

.026949 

.027474 

.028118 

0 028907 

.029872 

20 

.024390 

.024419 

.024507 

.024656 

. 024869 

.025150 

.025508 

.025951 

.026492 

.027149 

.0279^3 

21 

.023256 

.023281 

.023357 

.023486 

.023669 

„ 023912 

.024219 

.024597 

.025056 

.025609 

.026273 

22 

.022222 

.022244 

.022311 

.022423 

.022582 

.022792 

.023058 

.023383 

.023777 

.024248 

.024809 

23 

.021277 

.021296 

.021354 

.021452 

.021592 

.021778 

.022006 

.022289 

.022629 

.023034 

.023513 

24 

.020408 

.020425 

.020477 

.020563 

. 020686 

.020847 

.021049 

.021296 

.021592 

.021943 

.022356 

25 

.019608 

.019623 

.019668 

.019745 

.019854 

.019996 

.020175 

.020392 

.020651 

.020957 

.021316 

2  6 

.018868 

.018882 

.018922 

.018990 

.019087 

.019213 

.019371 

.019563 

.019792 

.020061 

.020375 

27 

.018182 

.018194 

.018230 

.018291 

.018377 

.018490 

.018631 

.018801 

.019004 

.019242 

.019519 

28 

.017544 

.017555 

.017587 

.017642 

.017719 

.017820 

.017946 

.018098 

.018279 

.018490 

.018736 

29 

.016949 

.016959 

.016988 

.017037 

.017107 

.017198 

.017311 

.017448 

.017609 

.017798 

.018016 

30 

.016393 

.016402 

.016429 

.016473 

.016536 

.016618 

.016720 

.016843 

.016989 

.017158 

.017353 

31 

.015873 

.015881 

.015905 

.015945 

.016003 

.016077 

.016169 

. 016280 

.016411 

.016564 

.016739 

32 

.015385 

.015392 

.015414 

.015451 

.015502 

.015570 

.015654 

.015754 

.015873 

.016011 

.016169 

33 

.014925 

.014932 

.014952 

.014986 

.015033 

.015094 

.015171 

.015262 

.015370 

.015495 

.015639 

34 

.014493 

.014499 

.014517 

.014546 

.014591 

.014648 

.014717 

.014801 

.014899 

.015013 

.015143 

35 

.014085 

.  014090 

.014107 

.014135 

.014175 

.014226 

.014290 

.014367 

.014456 

.014560 

.014679 

36 

.013699 

.013704 

.013719 

.013745 

.013782 

.013829 

.013888 

.013958 

.Ol4o4o 

.014135 

.014244 

37 

.013333 

.013338 

.013352 

.013376 

.013410 

.013453 

.013507 

.013572 

.013648 

.013735 

.013834 

38 

.012987 

.012991 

.013005 

.013027 

.013058 

.013098 

.013148 

.013207 

.013277 

.013357 

.013449 

39 

.012658 

.012662 

.012675 

.012695 

.012724 

.012761 

.012807 

.012862 

.012926 

.013000 

.013085 

40 

.012346 

.012349 

.012361 

.012380 

.012406 

.012441 

.012483 

.012534 

.012594 

.012662 

0 012740 

4l 

.0120*48 

.012052 

.012062 

.012080 

.012105 

.012137 

.012176 

.012223 

.012279 

.012342 

. 012414 

42 

.011765 

.011768 

.011778 

. 011794 

.011817 

.011847 

.011884 

.011928 

.011979 

.012038 

.012105 

43 

.011494 

.011497 

.011506 

.011522 

.011543 

.011571 

.011605 

.011646 

.011694 

.011749 

.011811 

44 

.011236 

.011239 

.011247 

.011262 

.011282 

.011308 

.011340 

.011378 

.011422 

.011473 

.011531 

45 

.010989 

.010992 

.011000 

.011013 

.011032 

.011056 

.011086 

.011121 

.OIII63 

.011210 

.011265 

46 

.010753 

.010755 

.010763 

.010775 

.010793 

.010815 

.010843 

.010877 

.010915 

.010960 

.011010 

47 

.010526 

.010529 

.010536 

.010547 

.010564 

.010585 

.010611 

.010643 

.010679 

.010721 

.010768 

48 

.010309 

.010311 

.010318 

.010329 

.010345 

.010365 

.010389 

.010418 

.010452 

.010492 

.010536 

49 

.010101 

.010103 

.010109 

.010120 

.010134 

.010153 

.010176 

.010204 

.010236 

.010272 

.010314 

50 

.009901 

.009901 

.009911 

.009921 

.009931 

.009950 

.009970 

.0099901 

.010030 

.OIOO69 

.010099 

k  / 
/  i 

0 

1 

2 

3 

4  * 

5 

6 

7 

8 

9 

10 

22 


MCA  TN  2173 


TABLE  I. -  COEFFICIENTS  \f±  -  Continued 


11 

12 

13 

14 

15 

1 6 

17 

18 
19 

0.20851 

.085295  0.20000 
.064150  .081893 
.053027  .061686 
.045836  .051065 

.040683  .044199 
.036752  .039281 
.033625  .035527 
.031061  .032541 

0.19245 

.078868 

.059486 

.049305 

.042727 

.038014 

.034417 

0.18570 

.076155  0.17961 

.057506  .073704  0.17408 

.047716  .055711  .071476  0.16903 
.041392  .046271  .054074  .069438 
.036862  .040176  .044950  .052574 

0.16440 

.067565  0.16013 

20 

.028909 

.030090 

.031555 

.033*105 

.035809 

.039060 

.043736 

.051192 

.065836  0.15617 

21 

.027070 

.028033 

.029205 

.030653 

.032478 

.034842 

.038033 

.042615 

.049913 

.064233 

22 

.025477 

.026274 

.027232 

.028394 

.029825 

.031623 

.033949 

.037082 

.041576 

.048726 

23 

.024079 

.024749 

.025544 

.026496 

.027647 

.029060 

.030833 

.033121 

.036200 

.o4o6io 

2k 

.022841 

.023410 

.024080 

.024872 

.025817 

.026956 

.028352 

.030099 

.032351 

.035378 

25 

.021735 

.022224 

.022794 

.023462 

.024250 

.025187 

.026314 

.027692 

.029415 

.031632 

26 

.020740 

.021163 

.021654 

.022224 

.022890 

.023673 

.024601 

.025716 

.027077 

.028776 

27 

.019839 

.020208 

.020634 

.021126 

.021695 

.022358 

.023135 

.024055 

.025157 

.026501 

28 

.019018 

.019343 

.019715 

.020143 

.020635 

.021202 

.021861 

.022632 

.023543 

.024633 

29 

.018267 

.018554 

.018882 

.019257 

.019685 

.020176 

.020741 

.021396 

.022161 

.023062 

30 

.017577 

.017832 

.018123 

.018453 

.018829 

.019257 

.019747 

.020309 

.020959 

.021717 

31 

.016940 

.017168 

.017427 

.017720 

.018052 

.018428 

.018856 

.019344 

.019903 

.020548 

32 

.016350 

.016555 

.016786 

.017048 

.017343 

.017676 

.018052 

.018479 

.018964 

.019520 

33 

.015802 

.015987 

.016195 

.016429 

.016693 

.OI6989 

.017322 

.017698 

.018123 

.018607 

3k 

.015291 

.015458 

.015646 

.015857 

.016094 

.016359 

.016656 

.016989 

.017364 

.017788 

35 

.014814 

.014966 

.015136 

.015327 

.015540 

.015778 

.016044 

.016342 

.016675 

.017049 

36 

.014367 

.014505 

.014660 

.014833 

.015026 

.015241 

.  015*181 

.015747 

.016045 

.016377 

37 

.013947 

.014073 

.014215 

.014373 

.014548 

.014743 

.014959 

.015199 

.015466 

.015763 

38 

.013552 

.013668 

.013798 

.013942 

.014102 

.014279 

.014475 

.014692 

.014932 

.015199 

39 

.013180 

.013286 

.013405 

.013537 

.013683 

.013845 

.014024 

.014221 

.014439 

.014679 

ko 

.012828 

.012926 

.013036 

.013157 

.013291 

.013439 

.013602 

.013782 

.013980 

.014198 

kl 

.012495 

.012586 

.012687 

.012799 

.012922 

.013058 

.013207 

.013372 

.013552 

.013751 

k2 

.012180 

.012264 

.012357 

.01^460 

.018574 

.012699 

.012837 

.012987 

.013152 

.013334 

4?1 

.011880 

.011958 

.012045 

.012140 

.012246 

.012361 

.012487 

.012626 

.012778 

.012944 

.011596 

.011668 

.011749 

.011837 

.011934 

.012041 

.012158 

.012286 

.012426 

.012578 

45 

.011325 

.011392 

.011467 

.011549 

.on64o 

.011739 

.011847 

.011965 

.012094 

.012235 

k6 

.OHO67 

.011130 

.011199 

.011276 

.011360 

.011452 

.011552 

.011662 

.011781 

.011911 

kl 

.010821 

.010879 

.010944 

.011016 

.011094 

.011180 

.011273 

.011375 

.011485 

.011605 

kQ 

.010585 

.010640 

.010701 

.010768 

.oio84i 

.010921 

.011008 

.011102 

.011205 

.011316 

lj-9 

.010360 

.  010412 

.010468 

.010531 

.010599 

.010674 

.010755 

.010843 

.010939 

.011043 

50 

.010149 

.010188 

.010248 

.010307 

.010366 

.010436 

.010515 

.010594 

.010693 

.010782 

,11 

12 

13 

l4 

15 

16 

17 

18 

19 

20 
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TABLE  I.  -  COEFFICIENTS  Akjl  -  Continued 


21 

0.15250 

22 

.062742  0.14907 

23 

.047620 

.061350  0.14586 

24 

.039708 

.046585 

.060047  0.14286 

25 

.034609 

.038864 

.045616 

.058824  0.14003 

*2  6 

.030959 

.033888 

.038071 

.044704  .057672  0.13736 

27 

.028176 

.030328 

.033211 

.037325  .043845 

.056585  0.13484 

28 

.025960 

.027613 

.029733 

.032572  .036621 

.043034 

.055558  0.13245 

29 

.024141 

.025452 

.027082 

.029173  .031969 

.035956 

.042266  .054585 

0.13019 

30 

.022610 

.023676 

.024971 

.026581  .028642 

.031399 

.035325  .041538 

.053661  0.12804 

31 

.021300 

.022184 

.023238 

.024517  .026106 

.028140 

.030857  .034727 

,o4o846 

.052782 

32 

.020161 

.020905 

.021780 

.022823  .024087 

.025656 

.027663  .030343 

.034158 

.040188 

33 

.019159 

.019794 

.020532 

.021398  .022430 

.023679 

.025229  .027210 

.029854 

.033616 

3^ 

.018268 

.018817 

.019447 

.020178  .021036 

.022056 

.023291  .024822 

.026778 

.029388 

35 

.017470 

.017946 

.018492 

.019117  .019842 

.020691 

.021700  *  .022921 

.024434 

.026366 

36 

.016749 

.017169 

.017644 

.018184  .018804 

.019522 

.020362  .021361 

.022569 

.024064 

37 

.016095 

. 016466 

.016883 

.017354  .017891 

.018506 

.019217  .020049 

.021037 

.022232 

38 

.015496 

.015826 

.016196 

.016611  .017079 

.017611 

.018221  .018926 

.019750 

.020728 

39 

.014946 

.015242 

.015571 

.015939  .016351 

.016816 

.017345  .017949 

.018647 

.019464 

ko 

.014439 

.014705 

.015000 

.015328  .015693 

.016103 

.016565  .017090 

.017689 

.018381 

kl 

.013969 

.014210 

.014475 

.014769  .015095 

.015459 

.015867  .016325 

-.016845 

.017440 

42 

.013532 

.013751 

.013991 

.014256  .014549 

.014873 

.015235  .015640 

.016095 

. 016611 

43 

.013125 

.013324 

.0135^3 

.013782  .014046 

.014338 

.Ol466l  .015021 

.015423 

.015875 

44 

.012745 

.012926 

.013125 

.013343  .013582 

.013846 

.014136  .014457 

.014815 

.015214 

45 

.012387 

.012554 

.012736 

.012935  .013153 

.013391 

.013653  .013942 

.014262 

.014617 

46 

.012052 

.012205 

.012372 

.012554  .012753 

.012970 

.013208  .013469 

.013756 

.014074 

47 

.011736 

. 011877 

.012031 

.012198  .012380 

.012578 

.012795  .013032 

.013292 

.013578 

48 

.011437 

.011568 

.011710 

.011864  .012031 

.012213 

.012411  .012626 

.012862 

.013121 

49 

.011155 

.011276 

.011407 

.011549  .011704 

.011871 

.012052  .012249 

.012464 

.012700 

50 

.OIO891 

.011000 

.011119 

.011257  .011396 

.011554 

.011723  .011891 

.012089 ' 

.012307 

k  X 
/  i 

21 

22 

23 

24  25 

26 

27  28 

29 

30 
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TABLE  I.-  COEFFICIENTS  Ak^  -  Continued 


31 

0.12599 

32 

.051945 

0.12404 

33 

.039560 

.051147 

0.12217 

3^ 

.033099 

.038961 

.050385 

0.12039 

35 

.028943 

.032605 

.038389 

.049655  0.11868 

36 

.025973 

.028517 

.032133 

.037841 

.048957  0.11704 

37 

.023710 

.025596 

.028110 

.031681 

.037315 

.048287  0.11547 

38 

.021910 

.023371 

.025236 

.027720 

.031247 

.036811 

.047644  0.11396 

39 

.020432 

.021601 

.023047 

.024890 

.027345 

.030830 

.036327 

.047026  0.11251 

40 

.019190 

.020148 

.021305 

.022735 

.024559 

.026985 

.030430 

.035862 

.046431  0.11111 

,  4l 

.018126 

.018927 

.019876 

.021021 

.022436 

.024240 

.026640 

.030045 

.035414 

.045858 

42 

.017201 

.017881 

.018675 

.019614 

.020748 

.022149 

.023933 

.026307 

.029674 

.034982 

43 

.016387 

. 016972 

.017646 

.018432 

.019363 

.020486 

.021872 

.023637 

.025986 

.029317 

44 

.015663 

.016171 

.016751 

.017419 

.018199 

.019121 

.020233 

.021605 

.023353 

.025677 

45 

.015014 

.015460 

.015964 

.016539 

.017202 

.017974 

.018888 

.019989 

.021348 

.023078 

46 

.014427 

.014822 

.015264 

.015764 

.016335 

.016992 

.017758 

.018663 

.019754 

.021100 

47 

.013894 

.014245 

.014636 

.015076 

.015572 

.016138 

.016790 

.017549 

.018446 

.019527 

46 

.013406 

.013720 

.014069 

.014458 

.014894 

.015387 

.015948 

.016594 

.017347 

.018236 

49 

.012957 

.013240 

.013553 

.013899 

.014286 

.014719 

.015208 

.015765 

.016406 

.017152 

50 

.012545 

.012792 

.013079 

.013386 

.013733 

.014119 

.014554 

.015035 

.015584 

.016224 

k  / 
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TABLE  I.-  COEFFICIENTS  ±  -  Concluded 


41  0.10976 

42  .045306  0.10847 

43  .034566  .044774  0.10721 

44  .028972  .034164  .044260  0.10600 

45  .025378  .028639  .033776  .043763  0.10483 

46  .022813  .025090  .028318  .033401  .043282  0,10370 

47  .020860  .022557  .024811  .028007  .033038  .042817  0.10260 

46  .019308  .020628  .022309  .024542  .027705  .032687  .042367  0.10153 

49  .018034  .019096  .020404  .022069  .024281  .027414  .032347  .041931  0.10050 

50  .016970  .017838  .018890  .020188  .021837  .024028  .027129  .032017  .041506  0.099504 
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(a)  Full  field  of  view 


Figure  1.-  Undisturbed  fringes  for  analyzing  interferogram  of  flow  at 

Mach  number  1.30. 
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(t>)  Region  used  in  analysis 


Figure  1.-  Concluded. 
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(a)  Full  field  of  view 


Figure  3--  Undisturbed  fringes  for  analyzing  interferogram  of  flow  at 

Mach  number  1.62. 
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Figure  3.-  Concluded. 
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(a)  Full  field,  of  view. 

V  L-6ij.901 

Figure  4.-  Interferogram  of  flow  around  sphere  at  Mach  number  1.62 


Figure  5.-  Contours  of  constant  ratio  of  density  to  free-stream  density. 


Figure  5.-  Concluded. 
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Figure  6.-  Variation  along  axis  of  ratio  p/pq- 
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Figure  7--  Variation  along  axis  of  ratio  p/ps. 


Local  Mach  number 


0  5  10  15  20  25  30  35  40  4  5  50  55  60  65  70  75  80  85  90 

Angular  distance  from  stagnation  point,  deg 

Figure  9*-  Variation  of  the  density  ratio  p/pq  with  angular  distance 
from  the  stagnation  point  along  the  surface  of  the  sphere.  M  =  1.6 
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Figure  10.-  Concentric  zones  in  disturbance. 
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Zone  number 


Figure  12.-  Assumed  distribution  of  density,  corresponding  distribution 
of  fringe  shift,  and  calculated  distributions  of  density. 
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NACA-Langley  -  8-23-50  -925 


Figure  13 Variation  of  ratio  p/pq  along  typical 


